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We report the discovery of energetic neutral hydrogen atoms emitted during the X9 solar 
event of December 5, 2006. Beginning ~1 hour following the onset of this E79 flare, the 
Low Energy Telescopes (LETs) on both the STEREO A and B spacecraft observed a 
sudden burst of 1.6 to 15 MeV protons beginning hours before the onset of the main solar 
energetic particle (SEP) event at Earth. More than 70% of these particles arrived from a 
longitude within ±10° of the Sun, consistent with the measurement resolution. The derived 
emission profile at the Sun had onset and peak times remarkably similar to the GOES soft 
X-ray profile and continued for more than an hour. The observed arrival directions and 
energy spectrum argue strongly that the particle events <5 MeV were due to energetic 
neutral hydrogen atoms (ENAs). To our knowledge, this is the first reported observation 
of ENA emission from a solar flare/coronal mass ejection. Possible origins for the 
production of ENAs in a large solar event are considered. We conclude that the observed 
ENAs were most likely produced in the high corona and that charge-transfer reactions 
between accelerated protons and partially-stripped coronal ions are an important source of 
ENAs in solar events. 

1 . Introduction 

In December of 2006, shortly after the launch of NASA’s twin STEREO 
spacecraft, there was a series of four X-class flares, each with an associated solar energetic 
particle (SEP) event. At this time the two STEREO spacecraft were located close together 
near Earth, and the Low Energy Telescopes (LETs) and High Energy Telescopes (HETs) 
were already operational. Measured SEP time-intensity profiles are shown in Figure 1. 
Since the December 5 event originated near the east limb, with no direct connection along 
the spiral interplanetary magnetic field, most energetic particles from this event arrived >4 
hours after the X-ray flare. Also seen between ~1130 - 1300 UT in Figure 1 is a small 
precursor to the main SEP event. This paper investigates the origin of this precursor. 

The 5 December event occurred with little warning and was not as well observed as 
most X-class events. There are no coronal mass ejection (CME) images because SOHO 
was undergoing routine maintenance and the STEREO/SECCHI instruments were not yet 
operational. However, as RHESSI emerged from Earth shadow at ~1030 UT it observed 
hard X-ray and y-ray emission near the peak in the GOES soft X-ray emission (see 
http://sprg.ssl.berkeley.edu/~tohban/browser/). In addition, at 1027 UT both STEREO 
WAVES instruments observed Type-II and Type-Ill radio bursts (showing that a shock 
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had formed in the corona and that electrons were escaping into interplanetary space; see 
http://swaves.gsfc.nasa.gov). 

2. Observations 

Together, LET and HET measure the nuclear charge (Z) and kinetic energy of H, 
He and heavier nuclei up through Ni from ~2 to ~100 MeV-nucleon. The LET detector 
system is composed of a double-ended array of 14 position-sensitive silicon solid-state 
detectors, including 10 thin (~25 micron thickness) detectors arranged in two fan-shaped 
arrays centered on a four-detector double-ended stack. The various detector combinations 
measure particle arrival directions over a 130° x 29° fan in both the front and back 
directions (Mewaldt et al. 2008). The HET sensor consists of a stack of ten circular silicon 
solid-state detectors (each 1 -mm thick) with a cone-shaped single-ended field of view with 
55° full angle (von Rosenvinge et al. 2008). 

The LET sensors on both spacecraft are normally oriented such that the Sunward fan 
is parallel to the ecliptic and centered 45° west of the Sun along the nominal (Parker spiral) 
interplanetary magnetic field. However, during December 2006 the STEREO-B spacecraft 
was flying up-side down such that the LET-B Sunward fan and HET-B were centered 45° 
east of the Sun. As a result, the two LETs combined could measure particle arrival 
directions over 360° of the ecliptic plane. 

The 5 December 2006 solar flare (located at E79 S06) was an X9 event, the sixth 
largest X-ray flare of solar cycle 23. The X-ray onset measured by NOAA’s GOES-12 
satellite occurred at 1019 UT with the maximum at -1030 UT. The arrival directions of 
1.6 to 15 MeV protons detected by the LETs during a 12-hour period on December 5 are 
shown in Figure 2. At ~1445 UT, protons, He, and heavier ions began arriving from a 
broad range of directions in the forward hemisphere. 

Also seen in Figure 2 at 1 130 - 1300 UT is a tight group of protons arriving over a 
very narrow angular range within ±10° of the Sun (no increase was seen in He or heavier- 
ion rates). The narrow width of the angular distribution and arrival times imply that these 
particles were not influenced by the magnetic field (oriented 36° to 67° from the Earth- 
Sun line during this interval) and we conclude that these protons must have been neutral 
particles (either neutrons or H atoms) while enroute to 1 AU that were converted to protons 
shortly before or as they were detected. 

Protons in this precursor exhibited velocity-dispersion with the highest-energy 
particles arriving first. Assuming these particles originated at the Sun and have not lost 
energy enroute, we derive the emission profile at the Sun using the measured kinetic 
energy (E) to determine the particle velocity, v = (2E/m) (using the proton mass for m). 
To compare to the X-ray time profile we plot the quantity T 0 b S - (Rst/v) + 492 seconds. 
Here Rst is the heliocentric distance of the STEREO spacecraft (0.983 AU). The added 
492 seconds correspond to the light travel-time from the Sun to the STEREOs (for X-rays, 
v is the speed of light and we plot the time profile reported by GOES). The resulting 
emission profile (Figure 3) is remarkably similar to the GOES X-ray profile, confirming 
that the neutral particles originated in this solar event. 

The distribution of the arrival directions is better illustrated in Figure 4, which 
includes 1.6-15 MeV protons arriving between 1015 and 1 145 UT. Note the very narrow 
peak in the arrival distribution within ±10° of the Sun. The typical maximum uncertainty 
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in these directions is ±6°, and most detector combinations that view the Sun were triggered 
during this burst. 


3. Discussion 

We first consider whether the LET observations could result from the decay of 
neutrons created when high-energy (>30 MeV) flare-accelerated protons and He 
underwent nuclear interactions with ambient H and He. This process has been modeled by 
Ramaty, Lingenfelter, and co-workers (Ramaty and Lingenfelter 1967, Hua and 
Lingenfelter 1987; Ramaty and Murphy 1987; Murphy et al. 2007). High-energy neutrons 
from a solar flare were first detected by Chupp et al. (1983, 1987). 

In addition, Evenson, Meyer and Pyle (1983; here-after EMP) and Evenson et al. 
(1990) measured interplanetary protons from the beta-decay of neutrons from three 
eastern-hemisphere events. Since the neutron mean-life is only 885 seconds, most neutrons 
emitted from the Sun decay well inside 1 AU. However, EMP showed that for east-limb 
events like this one neutrons that decay on field lines connected to Earth must have 
decayed relatively near Earth because of shadowing by the solar surface. 

In Figure 5a we compare the neutron-decay proton spectrum measured by EMP with 
energy spectra measured by LET and HET that include protons with derived solar emission 
times between 1015 and 1145 UT. All arrival directions are included and pre-event 
background proton rates were subtracted. Although neither HET viewed the Sun, there 
was a 2-o excess of 1 3 to 40 MeV protons during this time interval, and possibly a small 
excess of 40-100 MeV protons (shown as an upper limit). 

Also shown in Figure 5a is a theoretical neutron spectrum based on calculations by 
Hua and Lingenfelter (1987). Using their neutron spectrum emitted at 75° to the zenith we 
calculated the relative number of neutrons of each energy that would decay between 0.8 
and 1 AU and normalized the resulting decay spectrum to the LET and HET >5 MeV 
intensities. Although this estimate is approximate, it is clear that at 1 AU neutron-decay 
proton spectra of solar-flare origin should show a significant intensity decrease below ~10 
MeV due to the decay of low-energy neutrons before reaching field lines connected to 
STEREO. The LET spectrum below 5 MeV is inconsistent with an origin as neutron- 
decay protons. However, the excess 5 to 40 MeV protons following this flare may include 
a contribution from neutron-decay protons. 

While the initial direction of a neutron-decay proton will be aligned with the neutron 
trajectory from the Sun, the decay proton will immediately begin to spiral around the 
magnetic field. In addition, pitch-angle scattering will gradually isotropize the proton 
pitch-angle distribution (PAD; see Ruffolo, 1990), as observed by EMP. As a result, we 
would expect to detect decay protons from many directions. However, >70% of the <5 
MeV protons with emission times between 1015 and 1145 UT arrived from within 10° of 
the solar longitude, while at >5 MeV this fraction was ~23%. We conclude that neutrons 
were not responsible for the <5 MeV neutral emission, but could contribute at higher 
energy. 


Energetic Neutral Atoms (ENAs) are a second possible source of proton events in 
LET. Unlike neutron-decay protons, which originate at various points in the interplanetary 
medium both ahead and behind the spacecraft, ENAs originate at the Sun and travel in 
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straight lines. After leaving the solar vicinity most will remain neutral enroute to 1 AU. In 
LET they are stripped of their electrons upon striking the thin front Kapton windows. 
Solar ENAs thereby preserve their original directions until detection. 

ENAs have been observed from Earth’s magnetosphere and other planetary 
magnetospheres (Hovestadt and Scholer 1976; Decker et al. 2005; Mitchell et al. 2005) and 
also from the heliosphere (Hilchenbach et al. 1998). They are usually attributed to 
processes in which energetic protons (or heavier ions) suddenly become neutral through 
charge-exchange and then fly off in the same direction as that of the original charged 
particle. Most ENA observations have been made between ~1 and a few hundred keV. 

In Figure 5b we plot the fluence spectrum of the neutral burst, including only 
particles with arrival directions within ±10° of the Sun. Correcting the observed 1.8 to 5 
MeV ENA fluence back to the Sun, neglecting any losses and assuming isotropic emission, 
we estimate that 1 .8 x 10 ENAs escaped from the Sun in the upward hemisphere. 

The timing of the neutral-particle emission in Figure 3 is in excellent agreement 
with the X-ray flare, suggesting that the ENAs are produced by flare-accelerated particles. 
With ambient temperatures of 3 to 30 MK at the flare site (Feldman 1996), and 1 - 2 MK 
elsewhere in the corona, there should be no neutral H or He, and no He + , on which 
accelerated protons can charge exchange. However, neutral H atoms will also be produced 
by radiative recombination with ambient electrons (lF + e -> H + y) and by charge- 
exchange with heavier species such as O, Si, and Fe that retain some electrons at coronal 
temperatures. 

In the absence of measured or theoretical cross sections for charge-exchange reactions 
of protons with partially-stripped heavy ions at MeV energies, we have estimated the 
charge-transfer cross section for lF + 0 6+ -> H + 0 +7 from its inverse reaction, H + 0 +7 -> 
IF + 0 6+ , using the CTMC (classical trajectory Monte Carlo) method with detailed 
balancing (Olson and Salop 1977). There are some data for this inverse reaction <0.8 MeV 
(Phaneuf, Janev, and Pindzola 1987). While the simulation neglected corrections for K- 
shell capture by oxygen, which should enhance the high-energy yield (e.g., Stancil et al. 
1999), we consider this a reasonable first-order estimate for ENA production by proton 
collisions with O, the most abundant coronal species heavier than He. The CTMC result 
can be represented as cr = 1 .2 x 10‘ 20 E' 2 0 cm 2 from 0.1 to 5 MeV, with E the proton energy 
in MeV. We assume that cross sections for charge exchange with other coronal species 
scale from our oxygen estimate as (Z/8) 2 (Olson and Salop 1977). For radiative 
recombination we use On- = 1.28 x 10‘ 25 E' 20 cm 2 (with E in MeV) based on a fit to 
theoretical cross sections >1 MeV (Andersen and Boiko 1990), neglecting effects of 
electron heating in the flare region. 

Once produced, ENAs will be subject to ionization by electron and proton impact 
ionization and by UV. At MeV energies the electron and proton ionization cross sections 
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are equal and can be represented as Oj = 2.3 x 10’ x E' ' cm , with E the ENA energy in 
MeV (Barghouty 2000). 

In solar y-ray and neutron events ~10 31 to a few times 10 33 protons with >30 MeV 
are required to reproduce observations of y-rays and neutrons produced by interactions of 
flare-accelerated particles in the solar atmosphere (Ramaty and Murphy 1987). From 
Figure 5a we estimate that the 5 December flare may have produced ~10% as many 
neutrons as the 5 June 1982 event, for which Hua and Lingenfelter (1987) estimated that 
3.3 x 10 33 protons > 30 MeV were required. 
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Recent solar y-ray and neutron studies favor power-law spectra for flare-accelerated 
H and He with spectral indices ranging from -3.5 to -4.5 (Share and Murphy 1998; Lin et 
al. 2003; Murphy et al. 2007). For the 23 July 2002 event, Lin et al. (2003) derived proton 
spectral indices of -3.4 ±0.1, -3.8 ±0.1, and -4.0 ± 0.5 from three y-ray line ratios 
assuming a single power-law extending down to 2.5 MeV (the threshold for the Ne line). 

• 'ic • 

For a spectral index of -3.5 we estimate that ~3.5 x 10 protons in the range 1.8 to 5 MeV 
were accelerated in the 5 December 2006 flare event. 

We now estimate the ENA yield from this accelerated proton spectrum by 
assuming protons with <10 MeV slow down and stop somewhere in the solar atmosphere. 
Following Allison (1958), who analyzed the evolution of 2-component charge-exchange 
systems like H + and H, we calculate the mean number of >1.8-MeV ENAs created by a 
proton with energy E as it slows using: 

N mA = f*o c (E) [dE / dx(E)]~ 1 dE (1) 

Here Oc = OoiCfioAVoi + ciio), where aio is the abundance-weighted cross section for FT -» 
H (including radiative recombination and charge exchange), Ooi is the abundance-weighted 
cross section for the inverse reaction (electron and proton impact ionization), dE/dx(E) is 
the rate of energy loss in MeV/cm in an ionized plasma, and Nt is the number of target 
atoms per cm 3 assuming coronal abundances (Feldman and Widing 2003). Note that since 
Oio« Obi, Oc reduces to Oio. For a proton starting with 5 MeV we find an average yield of 
-0.01 ENAs before it slows to 1.8 MeV. In addition, >5 MeV protons can produce 1.8-5 
MeV ENAs after slowing to <5 MeV. Integrating over the adopted solar-flare spectrum 
from 1.8 to 10 MeV we get ~10 33 ENAs produced with 1.8-5 MeV. Assuming isotropic 
emission, this is ~30,000 times the number needed to explain our observations if all 
upward-moving ENAs escape the Sun (attenuation in the solar atmosphere is considered 
below). 

In the standard picture of a solar flare a great deal of energy is suddenly released in 
the corona through magnetic reconnection and energetic particles are quickly accelerated 
by one or more processes (e.g., Miller 1998; Aschwanden 2002). Based on studies of X- 
ray flares, Aschwanden (2002) concluded that the mean height of the acceleration region 
varies from ~5000 to 50,000 km, with an extreme maximum of ~200,000 km, 
corresponding to a heliocentric radius of ~1.3 Rs- The acceleration region produces both 
upward and downward-directed beams of energetic electrons and protons. Downward- 
directed electrons produce bremsstrahlung emission from the chromospheric footpoints of 
the flare, while nuclear interactions of downward-directed ions are responsible for solar y- 
ray and neutron events (e.g., Lin et al. 2003). 

Acceleration sites <1.3 Rs will be copious sources of ENA production, but when 
the overlying solar atmosphere in the direction of Earth is considered the estimated 
attenuation factor for a typical 3-MeV ENA produced at 1.15 to 1.3 Rs ranges from ~10" 24 
to 10’ . The attenuation factor is given by F(R) = exp(-OiNR) where Nr (in cm' ) is the 
overlying column density of protons and electrons integrated from heliocentric distance R 
(at E79) to the STEREOs, using nominal coronal densities (Sittler and Guhathakurta 
1999). This indicates that ENAs from the acceleration site or the chromospheric footpoints 
should not be detectable at Earth for east-limb flares. Indeed, for production sites <1.6 Rs 
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we find that not enough ENAs escape towards Earth in this event, even if all 10 ENAs 
were produced there. Higher in the corona the estimated 3-MeV ENA attenuation toward 
Earth ranges from ~10' 4 at 1 .7 Rs to ~10' 2 at 2.0 Rs. We suggest that ENA emission in the 
December 5 event resulted mostly from ENAs produced beyond ~1 .6 Rs. 

ENA production will also occur by particles accelerated at the CME-driven shock. 
Using the same cross sections, coronal densities, an E' 1 spectrum as observed at 1 AU 
(Mewaldt et al. 2005), and assuming all accelerated particles are concentrated near the 
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shock, we estimate that ~10 protons with 1.8-5 MeV are needed to account for our 
measured ENA fluence (Figs 3 and 5b) as the shock moves from 3 Rs to 9 Rs in 35 min 
(assuming a 2000 km/s shock). This number of protons is comparable to the estimated 
source strength of 1.8 - 5 MeV protons from the large 28 October 2003 event (Mewaldt et 
al. 2005), indicating that the ENA yield from CME shock-accelerated particles alone may 
be inadequate. However, more realistic estimates would include the lateral spread of the 
shock and ENA production in higher-density material behind the shock. 

These ENA yields depend on our estimated cross-sections for charge exchange 
with partially-stripped heavy ions, which are only first-order estimates. There are also 
uncertainties in the assumed solar-flare proton spectrum and in coronal densities. 

4. Summary and Conclusions 

This paper reports measurements of an energetic neutral-hydrogen burst from 
longitudes within ±10° of the Sun during the large 5 December 2006 solar flare. The 
derived solar emission profile is in excellent agreement with the soft X-ray profile if the 
particles traveled in a straight line from the Sun. The arrival directions, energy spectra, 
and timing of this neutral burst argue strongly that most of the particles (especially those 
<5 MeV) were energetic neutral hydrogen rather than neutrons. To our knowledge this is 
the first reported observation of ENAs >10 keV/nuc from a solar source. 

In a flare or coronal environment energetic protons produce ENAs by radiative 
recombination. In addition, our rough estimates suggest that charge-exchange reactions 
with partially-stripped heavy ions may produce ~100 times as many ENAs as radiative 
recombination. We estimate an upper limit to 1.8-5 MeV ENA production in this event 
of ~10 33 ENAs by assuming flare-accelerated protons <10 MeV slow and stop in the solar 
atmosphere. By comparison, our observations indicate that ~1.8 x 10 ENAs escaped the 
Sun (assuming isotropic emission). Taking into account ENA losses we conclude that 
ENAs from the reconnection site or flare foot-points would not be observable at Earth. 
Rather, observations of east-limb events heavily favor ENA production beyond ~1.6 Rs, 
including flare-accelerated protons that reach the high corona and particles accelerated at 
the CME-driven shock. More detailed models of ENA production and loss in a solar- 
flare/CME environment will be required to assess these and other possible ENA 
production sites. 

In the approaching solar maximum, STEREO ENA observations from multiple 
points of view, coupled with models and other imaging data, can provide a new window 
for studying particle acceleration and transport on the Sun by revealing where and when 
low-energy ions interact with solar material. In particular, solar ENA observations may 
provide a means of investigating the poorly-known low-energy spectra of ions accelerated 
in solar flares. 
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Day of December, 2006 



Fig. l.-Time history of low-energy protons measured by LET-B during 5 to 8 
December 2006. The onset time and location of two X-class flares are indicated. 
The arrow points to a burst of particles preceding the main SEP event. 
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Hour of 5 Dec 2006 


Fig. 2 -Plot of the measured angle to the Sun for individual 1.6 to 12 MeV proton 
events on 5 December 2006 (red = LET-A; blue = LET-B). Note the group of events 
within ±10° of the Sun from ~1130 UT to ~1350 UT, well before the SEP onset at 
~1445 UT. The X9 flare onset is indicated. The range of magnetic field orientations 
connecting to the Sun between 1130 and 1350 UT is shown as a bar. 
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Fig. 3 -The derived emission profile of the 1.6 - 15 MeV neutral particle burst (in 
counts per 5 minutes) is compared with the GOES X-ray profile and with the onset of 
STEREO type-II and type-III radio bursts. The neutral-particle and X-ray emission 
peaks agree within ±2.5 minutes. Both traces have 5-minute resolution. 
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Arrival Direction (Sun = 0°) 


Fig. 4.-Histogram of arrival directions for LET proton events that arrived during 
1015 - 1145 UT on December 5. The angular coverage of LET- A (solid) and LET-B 
(dashed) are indicated along the top. 
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Fig. 5. -(a) Proton spectra measured by LET and HET (including protons from all 
arrival directions if the derived solar emission times are between 1015 and 1 145 UT 
on December 5) are compared to neutron-decay proton spectra from IMP-8 for the 3 
June 1982 event (Evenson et al. 1983). Also shown (dotted line) is a 1-AU neutron 
spectrum based on calculations by Hua and Lingenfelter (1987; see text). Five LET 
events with 1.6 to 1.8 MeV were excluded because of detection-efficiency 
uncertainties, (b) Spectrum of energetic neutral hydrogen measured by LET- A and 
LET-B during the 5 December 2006 solar event, based on detected protons with 
arrival directions within ±10° of the Sun and derived solar emission times between 
1015 and 1145 UT. The >5 MeV points may also include a contribution from 
neutron-decay protons (see text). Note the different units in Fig. 5a and Fig. 5b. 


- 12 - 


